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1. Introduction 

In the last few years remarkable mathematical structures were shown to emerge in the 
analysis of supersymmetric Chern-Simons-matter theories in three dimensions. In this 
class of theories a distinguished role is played by the Af = 6 ABJM model [1] which is 
a U(N) k x U(N)_ k superconformal gauge theory with Chern-Simons level k. Indeed, 
in the large N limit the ABJM theory has been conjectured to be the AdS/CFT dual 
description of M-theory on an AdS± x Sj/Zk background and, for k <C iV <C k 5 , of 
a type IIA string theory on AdS^ x CP3. For this reason, soon after its discovery 
the ABJM model has quickly become the ideal three-dimensional playground to study 
AdS/CFT as much as Af = 4 SYM has been in the four-dimensional case. 

Quite surprisingly, the ABJM model seems to share a number of remarkable prop- 
erties with Af = 4 SYM theory even if the two theories are a priori different in nature. 
One of the common features is provided by the fact that also in the ABJM case, in the 
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limit of large N and k with A = N/ k kept fixed, integrable structures naturally show 
up (for a review see [2]). 

At first it was found in [3, 4] that, at the two-loop order and in the SU(4) fla- 
vor sector, the anomalous dimensions of composite operators could be mapped to the 
energy spectrum of an integrable Hamiltonian acting on an alternating fundamental- 
ant if undamental spin-chain. 

The two-loop analysis was then extended to the full theory in [5, 6] by the intro- 
duction of an OSp(2, 2|6) chain. The parity breaking ABJ model [7] was also studied 
at two loops in [6, 8], where it was found to be integrable at the given order. 

Soon afterwards, paralleling the progresses done in the four-dimensional set 
of all-loop Bethe equations for the asymptotic spectrum of the full ABJM theory was 
proposed in [9]. The Bethe equations nicely interpolated between the weak coupling 
results and the coset string construction at strong coupling [10, 11], together with the 
algebraic curve approach developed in [12]. 

One of the salient features of the Bethe equations introduced by Gromov and 
Vieira is that they are strongly constrained by the symmetries of the theory. In fact 
only a pair of undetermined functions of the coupling A are left open in the description 
of the spectrum. One is the interpolating function h(X) which relates the weak and 
strong coupling regimes of the single magnon dispersion relation. The other one is 
the dressing function 0(A), which had to be introduced also in the M = 4 SYM case 
[13, 14, 15, 16, 17, 18] in order to have a proper matching between the weak and 
strong coupling descriptions. In the ABJM case, the dressing phase plays an even more 
fundamental role since its presence has also been conjectured in [9] to give rise to the 
coupling between even- and odd-site excitations on the chain at high-loop orders. 

The description of the asymptotic spectrum depicted in [9] has been subsequently 
checked at the perturbative level beyond two-loops by direct Feynman diagrammatic 
computations in [19, 20, 21, 22]. The four-loop dilatation operator has been fully 
computed for both the ABJM and ABJ models by using the component formulation 
in [19, 20] and the M = 2 superspace formalism in [21]. As a result, the form of 
dilatation operator was found to be compatible with the spectrum predicted by the 
Bethe equations and moreover it was possible to fix the next-to-leading order coefficient 
in the weak coupling expansion of the function h(X). 

In [23] the analysis was pushed up to six loops. At this order, a full Feynman 
diagrammatic analysis looks very complicated, even using superspace techniques. Nev- 
ertheless, in [23] the expression of the dilatation operator could be derived by computing 
a suitable set of Feynman diagrams. For further perturbative checks on the spectrum 
see [24]. 

Meanwhile, the internal S-matrix approach to integrability was developed in a 
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number of papers [25, 26, 27, 28, 29, 30]. An all loop S-matrix has been found which 
is compatible with the all-loop Bethe equations of [9]. Moreover, the thermodinamic 
Bethe Ansatz and Y-system framework has also been applied to the three dimensional 
case [31, 32, 33, 34]. 

In the present paper we would like to provide further evidence for the exactness 
of the conjectured integrable scenario by performing a direct computation of the first 
non- vanishing coefficient of the Bethe Ansatz dressing phase. We restrict ourselves to 
the SU (2) x SU (2) scalar sector of the theory, where the spin chain is simply given by 
a pair of coupled SU (2) spin chains. In [9] it was mentioned that the dressing phase 
starts contributing at eight loops, leading to the result that the two SU (2) chains at odd 
and even sites are decoupled up to six loops. Indeed, it was explicitly shown in [19, 20] 
and [21] by direct computations, that this is true at four-loop order: the contributions 
to the dilatation operator of the diagrams that could lead to interactions between the 
two types of magnons cancel. Moreover, the results of [23] imply that this is true also 
at six-loop order. This latter fact isn't so trivial and has the consequence that the 
order A 4 coefficient of the dressing factor which can be a priori present at six loops, 
actually vanishes. Therefore, the first non-trivial effects of the presence of the dressing 
phase are to be found at eight-loop order. With the aim to check the above picture, 
we analyze the form of the eight-loop dilatation operator and extract the value of the 
leading order coefficient of the dressing factor. We make use of a procedure inspired to 
the one used in [35] for the computation of the analogous coefficient in J\f = 4 SYM. 

The plan of the paper is the following. After reviewing in Section 2 the aspects of 
integrability in the ABJM model we shall need, in Section 3 we introduce a procedure 
that is useful to constrain the form of the dilatation operator as much as possible. 
This procedure is based on symmetry arguments and on a matching between the en- 
tries of the diagonalized dilatation matrix with the Bethe equations predictions for 
the energies of one and two-impurity states. We then use this procedure to write the 
maximal reshuffling part of the dilatation operator in terms of a small number of un- 
known parameters, one of which is the dressing phase. In Section 4, we perform a 
direct Feynman diagram computation of the maximal reshuffling diagrams to fix the 
unknown coefficients including the dressing phase. We end up with a result which is in 
complete agreement with the integrable picture conjectured in [9]. As we shall discuss, 
our computation, besides providing the explicit value for the leading coefficient of the 
dressing factor, represents a non-trivial consistency check of the Bethe equations, since, 
as mentioned before, the dressing phase is also shown to play the structural role of cou- 
pling even and odd excitations in the ABJM model. We give further comments on our 
results in Section 5 while several technical aspects are collected in the Appendices. 
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2. The Bethe Ansatz in the 577(2) x 577(2) Sector 

In this section, we collect the integrability tools we need for the construction of the 
dilatation operator. We refer to Appendix A for the formulation of ABJM theory in 
terms of M = 2 superfields. Throughout the entire work, we restrict ourselves to the 
577(2) x SU(2) sector, where operators are made out only of the chiral superfields Z A , 
W B : 

Oi\::£ = Tr[Z Al W Bl ■ ■ ■ Z a ^W Bl ] , (2.1) 

with Aj, Bj = 1, 2. This sector is closed under renormalization at all loop-order. In the 
spin-chain picture, the operators (2.1) are mapped to states of a circular alternating 
spin-chain, with the fields Z A being interpreted as spins lying on odd sites and the Wa 
as spins lying on even sites. We will refer to fields Z 2 and W 2 as impurities and they 
correspond to "spin down" states. The ground state of length 2L is chosen to be 

| ) = Tr^Wi ■ ■ ■ Z X W^ , (2.2) 

while the operators 

\2k) =Tr[Z 1 W 1 ---Z 1 W 2 ---Z 1 W 1 ], 
\2k + 1 ) = TrfZ 1 ^ • • • Z 2 W X ■ ■ ■ Z X W X \ 

represent states with a single excitation on the site 2k and 2k + 1 respectively. There 
are two kinds of magnon states, depending on their momentum p being excited on even 
or odd sites 1 : 



|p) e = 5>**|2*> , 

k=l 
L 

\p) = ^e ipk \2k + l) 



k t (2-4) 



k=i 

The magnon dispersion relation is given by 



E{p)=\ ^l + 16^(A)sin 2 |-lj , (2.5) 
where h(X) is the interpolating function, which has the weak-coupling expansion 2 

oo 

h\\) = h 2k \ 2k (2.6) 
k=i 



1 In the sum, the identification 2L + 1 ~ 1 is understood. 

2 For discussions on the strong coupling expansion of h(X) see [36]- [45]. 
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with Y12 — 1 and h± = — 4£(2) the only known coefficients [19, 20, 21]. 

The all-loop asymptotic Bethe equations for the length 2L spin chain can be ob- 
tained from [9] by restricting to the SU(2) x SU (2) flavor sector 



x(uj + i/2) 
x(uj — 

' x(vj+i/2) 
x(vj - i/2) 



l L 



M u 

n 



k=1M «i ~u k -i 



l L 



M v 

n 



M v 

n 

k=i 



,id(uj,v k ) 



.Vj-Vk-i 



k=l,k^j J k=l 

and must be supplemented by the momentum constraint 



yT xjUj+i/2) yr xjvj+i/2) 
f = \x(v 3 -t/2) 



1. 



(2.7) 



(2.8) 



We have called Uj and Vi the Bethe roots of each SU(2) factor. We also denoted by M u 
and M v the number of u and v roots. Moreover, we introduced the function 



w 



x(w) = -[ l + 



h?{\) 



(2.9) 



where w = u,v. As we see from (2.7), the dressing factor e ld introduces extra self- 
interactions for the roots and also couples the two SU(2) spin chains at higher loops. 
Following [9] we define the dressing phase as in [14, 17] 



oo oo 



r=2 s=r+l 

where the coefficient functions /? rjS (A) can be expanded in A as 



(2.10) 



MA) = ti 2 s )x2k 



k=s-l 



and the quantities 



q r (w) 



r - 1 \x(w + i/2) 7 - 1 x(w - i/2) r - 1 



(2.11) 



(2.12) 



are related to the eigenvalues of the conserved charges of the theory, whose existence 
is ensured by integrability. In particular, the second charge is the hamiltonian of the 
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integrable model and is identified with the dilatation operator of the gauge theory. The 
energy eigenvalues are given by 



+ • (2-13) 

3=1 j=l J 

Let's observe that, for two-impurity states, the Bethe Ansatz is simplified. The mo- 
mentum constraint requires 

w 1 = —w 2 = w , (2-14) 
while the Bethe equations reduce to 



w + i/2 \ L _ ( 1 + V 1 ~ FW \ pi6{w ,- w) 



W—i/2) \l_i_ ft 4fe 2 (A) 

\ 1 + V (i«+i/2) 2 



e iv{w-w } (2.15) 



for M u = M v = 1, and 



w + i/2\ 1 '" 1 / 1 + \/ 1_ (^=W 



L 

= j0(io,-'iu) 



W — i/2 J \1i /-i 4fe 2 (A~ 

V 1 + V (io+i/2) 2 



for M u = 2, M v = or M u = 0, M„ = 2. For any fixed L, these equations can be solved 
order by order in perturbation theory. This is the basic tool for the computation of 
anomalous dimensions of long operators. 

As can be seen from the Bethe equations, the anomalous dimensions of operators 
will depend, in general, on the coefficients of the function h(X) and of the dressing 
phase. It seems that such coefficients cannot be determined using the integrability of 
the theory only. However, some general considerations on their values can be made, as 
we now explain. In the case of M = 4 SYM theory it was shown [17] that the coefficients 
of the dressing phase should have a well defined degree of transcendentality in order to 
preserve the Kotikov-Lipatov transcendentality principle on the scaling function of the 
theory [46]. Such principle can be generalized to the ABJM case: the scaling function 
of ABJM theory is prescribed to be [9] 

/abjm(A) = - /sym(sO| s -^(a) , (2-17) 
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where g = ^ in M = 4 SYM. Taking the weak coupling expansion for /sym(<?) obtained 
in [17] and inserting the expansion of h(X), we thus get, up to eight loops, 



/abjm(A)=4A 2 




V - 4fc 4 ) A 4 + ^ - 8 -h^ + 4/*) A 6 

- 8C(3) 2 - ^ 4 vr 4 + \hW + -h 8 + 4/3gC(3)) A' 



8 



(2.18) 



In deriving (2.18) we have used the fact that (3 TyS = for r + s even and that, in 
our case, the coefficients /3 rjS (A) should contribute starting at order 0(A 2(r+s ~ 2 )) rather 
than (9(A 2(s_1 )), as expected from (2.11). See [14, 17] for an explanation. The latter 
property will be verified, for the /3 2 ,3 coefficient, in Section 3. The transcendentality 
principle, extended to ABJM theory, states that: 

Assigning degree of transcendentality k to constants n k and((k), the l-loop contribution 
to the scaling function f abjm has uniform degree of transcendentality t — 2. 

Recalling that h 4 = —4^(2), we see that this principle is satisfied up to four loops. 
Moreover, from the six-loop contribution, we immediately see that he should have de- 
gree four. Starting from eight loops also the unknown dressing phase starts contributing 
to the scaling function. One very natural way 3 to preserve the transcendentality prin- 
ciple is to conjecture that: 

• h 2 k has degree of transcendentality 2k — 2, 

• has degree of transcendentality 2k + 2 — r — s. 

In particular, we notice that the transcendentality principle implies that the leading 
order coefficient of the dressing phase should be of the form 



where a and b are some rational numbers. The main goal of this paper is to compute 
such constants. 




avr 3 + K(3), 



(2.19) 



3 We're neglecting the possibility that hs and p^. 3 can have higher (or lower) degrees of transcen- 
dentality which cancel out between them. 
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3. Construction of the Dilatation Operator 



The dilatation operator is the generator of scaling transformations and measures the 
scaling dimensions of composite operators. The perturbative expansion of the dilatation 
operator in the SU (2) x SU (2) sector is 

oo 

P(A) = L + ^P 2fc A 2fc , (3.1) 
k=\ 

where the 0-loop term L is proportional to the identity and yields the classical dimension 
of an operator made up of 2L chiral superfields, while the second term is the quantum 
contribution and yields its anomalous dimension. The latter can be extracted from the 
perturbative renormalization of the composite operators O a 

£>r = Z a b O h b we , Z = t + \ 2 Z 2 + \ 4 Z 4 + ■■■ . (3.2) 

The renormalization factor Z is introduced in order to remove UV divergences from 
correlation functions of operators and can be computed in perturbation theory by 
means of standard methods. Using dimensional regularization in spacetime dimension 
D = 3 — 2e, quantum divergences show up as inverse powers of e in the limit e — > 0. By 
introducing the 't Hooft mass \x and the dimensionful combination A/z 2e the dilatation 
operator is then extracted from Z as 



d 

V = L + a— In Z(\a 2£ , e) = L + lim 

d/J, e^O 



2e\^-\nZ(\,e) 
clA 



(3.3) 



This definition effectively extracts the 1/e pole of InZ; the higher order poles must be 
absent. A full-fledged quantum field theory computation of the dilatation operator can 
be an extremely difficult task, especially for high loop orders in perturbation theory. 
Fortunately, things get simplified if we make use of the integrability results. 

Here we discuss a general procedure to construct the asymptotic dilatation operator 
of ABJM theory. It is similar to the procedure used in the case of A/" = 4 SYM 
and described in [47, 48, 49, 14, 50]. The idea is the following: starting from an 
exactly integrable hamiltonian (which is next-to-nearest neighbor for ABJM theory), 
we deform it with local interactions of range linearly increasing with the perturbative 
order, obtaining a long-ranged spin chain. At I loops the maximum range is I + 1. 
When the range of an interaction exceeds the length 2L of the spin chain, wrapping 
interactions typically appear and the Bethe Ansatz breaks down. So, we work with 
asymptotic states, i.e. states of length 2L > £. 

The two-loop dilatation operator [3, 4] is the sum of two XXXi Heisenberg hamil- 
tonians, one living on odd and the other on even sites: this hamiltonian is deformed 
with interactions of longer range at higher perturbative orders. 
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The dilatation operator is conveniently written in terms of the basis of permutation 
structures 

L 

{ai, . . . , a n } = ^2 P2j+ ai ,2j+ ai +2 ■ ■ ■ P2j+a n ,2j+a n +2 , (3.4) 
J'=l 

where Pjj are permutation operators. Permutation structures represent local interac- 
tions of spins summed over all different positions on the spin chain: their properties 
are collected in Appendix B. Here, we just need to define the number 

1Z = max(ai, . . . , a n ) — min(ai, . . . , a n ) + 3 (3.5) 

which is the range of the interaction. 

In Table 1 we list all the independent permutation structures that appear in the 
dilatation operator up to eight loops. In order to read a basis for the dilatation operator 
at t loops, one has to include all the permutation structures up to range 1Z — £ + 1, 
neglecting the ones in parenthesis, because they appear at the next loop order. 

In order to construct the £-loop dilatation operator, we first write the most general 
combination of permutation structures up to range £ + 1. Then, we fix a large part of 
the unknown coefficients as follows: 

1. Impose hermiticity and parity invariance: 

v\ = v e , w/p- 1 = V t ■ 

2. Require that the vacuum state is protected (i.e., it has zero energy): 

V e \0) =0; 

3. Impose the magnon dispersion relation on one-impurity states on odd and even 
sites (let E e be the £-loop coefficient of (2.5) in the small A expansion): 

V t \p) 0:e = E e (p) \p) , e \ 

4. Use the asymptotic Bethe ansatz on two-impurity states in order to fix the re- 
maining parameters. 

Point 4 works in the following manner: 

• Fix a length of the spin chain, with L > £/2, so that wrapping interactions are 
absent; 
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TZ 


Basis of permutation structures 


1 


{} 


2 




3 


{0} 


4 


{0,1} 


5 


{0,2}, {2,0} 
({0,1,2}, {2,1,0}) 


6 


{0,3} 
{0,1,3}, {0,3,1} 
{0,2,3}, {2,0,3} 
({0,1,2,3}, {2,1,0,3}, {0,3,2,1}, {2,3,0,1}) 


7 


{0,4} 

{0,2,4}, {4,2,0}, {0,4,2}, {2,0,4} 
({0,1,4}, {0,3,4}) 
({2,0,4,2}) 

({0,1,2,4}, {2,1,0,4}, {4,1,2,0}, {4,1,0,2}) 

({0,1,3,4}, {0,3,1,4}) 
({0,2,3,4}, {2,0,3,4}, {0,4,3,2}, {4,3,2,0}) 


8 


{0,5}, {0,1,5}, {0,4,5} 
{0,2,5}, {2,0,5}, {0,3,5}, {0,5,3} 
{0,1,3,5}, {0,3,1,5}, {0,1,5,3}, {5,3,1,0} 
{0,2,3,5}, {2,0,3,5}, {0,2,5,3}, {2,0,5,3} 
{0,2,4,5}, {2,0,4,5}, {0,4,2,5}, {4,2,0,5} 


9 


{0,6} 

{0,2,6}, {2,0,6}, {0,4,6}, {0,6,4} 
{0,3,6} 

{0,2,4,6}, {0,2,6,4}, {0,4,2,6}, {2,0,4,6} 
{2,0,6,4}, {0,6,4,2}, {4,2,0,6}, {6,4,2,0} 



Table 1: Permutation structures needed up to eight loops, grouped according to their range 
TZ. Permutation structures of odd (even) range TZ appear at TZ — 1 (TZ) loops, apart from the 
ones in parenthesis, which appear at the next loop order. Permutation structures contributing 
at ten loops or beyond are not written. Note that, in order to find the complete basis, one 
has to add {. . . ,a + 1, . . .} to each permutation structure {. . . , a, . . .} listed here. 



• Find a basis of two-impurity operators of length 2L. Note that there are two 
types of such bases, as the impurities can be both on even (or odd) sites or one 
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on even and the other on odd sites 4 ; 

• Explicitly diagonalize the matrix representation of T>(\) 
basis of length 2L states; 

• Compute the eigenvalues from the asymptotic Bethe equations (2.7) and compare 
them with those found in the previous step, containing the unknown parameters. 

The explicit diagonalization of the dilatation operator and the solution of the Behte 
equations can be performed with the help of Mathematica. 

At the end of the procedure, some parameters may remain unfixed. Some of them 
are related to similarity transformations of the dilatation operator: 

V = e~ ix Ve ix . (3.6) 

They don't appear in the spectrum but only affect the eigenvectors. Therefore, they 
can be regarded as unphysical and their values depend on the renormalization scheme. 

Another way to fix the unknown parameters is through conservation of the first 
higher charge Q3, if we assume perturbative integrability. This requires the construction 
of Q 3 up to order £ with a procedure similar to that described for the dilatation operator, 
with the difference that has odd parity and is anti-hermitian. Since Q3 has higher 
maximal range, £ + 2, at high loops it can be quite complicated to compute. Some 
other parameters may depend on the dressing phase and cannot be fixed through the 
Bethe Ansatz. 

We now apply the procedure described above to constrain the form of the dilatation 
operator up to eigth loops. We immediately get 

V 2 = 2{}-{0}-{l} (3.7) 

for the two-loop dilatation operator and 

X? 4 = 2 (-4 + h 4 ) { } + (6 - h A ) ({0} + {1}) - ({0, 2} + {2, 0} + {1, 3} + {3, 1}) (3.8) 

for the four-loop dilatation operator. We see that (3.7) reproduces the result of [3, 4], 
while (3.8) reproduces the results 5 of [19] and [22]. 

4 For the same reason, we have seen in section 2 that there are two types of Bethe equations for 
two-impurity states. 

5 The computation of [19] followed a procedure similar to ours, up to point 3. This fixes all but one 
coefficient in the dilatation operator, namely the coefficient of {0, 1} + {1, 2}. We simply use point 4 
of our procedure on length six two-impurity states to prove that it vanishes, without computing any 
Feynman diagram. 



= Ek=i V 2k \ 2h on the 
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The six-loop dilatation operator was computed in [23] through Feynman diagram- 
matics. We apply our technique to rederive it. Assuming that the dressing phase is 
absent at this loop order, following [14, 17], we can completely fix the dilatation oper- 
ator with our procedure. If we don't assume that the dressing phase is absent, a single 
input from the results of [23] is sufficient to determine the full dilatation operator. 
We follow this latter approach. First of all we write the six-loop dilatation operator 
in the basis of permutation structures. After having imposed hermiticity and parity 
invariance, it is given by 

,mixed j 

(3.9) 

where 

V 6 , even = b {0} + d ({0, 2} + {2, 0}) + m {0, 4} 

+ I ({0, 2, 4} + {4, 2, 0}) + [({2, 0, 4} + {0, 4, 2}) (3.10) 
+ ^({2,0,4} -{0,4, 2}), 

T>^odd = ^6,even({- .., a, ...}**{..., a + 1, ...» , (3.11) 



V 6>mixed = c ({0, 1} + {1, 2}) + e ({0, 3} + {1, 4}) 

+ / ({0, 1, 2} + {2, 1, 0} + {1, 2, 3} + {3, 2, 1}) 
+ g ({0, 1,3} + {0, 3, 1} + {1, 3, 4} + {3, 1, 4} 

+ {0, 2, 3} + {2, 0, 3} + {1, 4, 2} + {1, 2, 4}) 
+ i e 2 ({0, 1, 3} - {0, 3, 1} - {1, 3, 4} + {3, 1, 4} 

- {0, 2, 3} + {2, 0, 3} - {1, 4, 2} + {1, 2, 4}) , 

where all parameters are real. Zero-energy of the vacuum gives the equation 

a + 2b + 4d + 41 + 41 + 2m + 2c + 2e + 4/ + 8# = . (3.13) 

We then compute the action of T>§ on \p) e and \p) and compare the result with the six 
loop coefficient of the weak-coupling expansion of (2.5) obtaining the equations 



/= -2, 

b = - 48 - 2c - 2e - 2f - 4g + I2h 4 - h 6 + 21- 2m , 
d = 12- f -2g-2h 4 -2l. 



(3.14) 



Now we consider two-impurity states. In order to avoid wrapping, we have to start 
from L — 4. For L = 4 (length 8 states), comparing the eigenvalues of the dilatation 
operator with the solutions of the Bethe equations we get 

[=0, 2m = 4- 3/9$, f = 0, g = 0, c + e = 0. (3.15) 

The results of [23] (obtained with an independent approach) imply, in our notations, 
that m = 2. From (3.15) we thus see that (3^1 must vanish. For L = 5 (length 10 

states) we get 

c = 0, e = 0. (3.16) 

We therefore have explicitly shown that, up to irrelevant parameters, T> &mixed = 0. It 
means that the two types of magnons don't interact at six loops, i.e. the two 577(2) 
factors are decoupled, as claimed in [9]. 

The six-loop dilatation operator is then given by 

T? 6 = 2 (30 - 8h 4 + h 6 ) { } + (-52 + 12h 4 - h 6 ) ({0} + {1}) 

+ 2 ({0, 4} + {1, 5}) + (12 - 2h A ) ({0, 2} + {2, 0} + {1, 3} + {3, 1}) 
- 2 ({0, 2, 4} + {4, 2, 0} + {1, 3, 5} + {5, 3, 1}) 
+ 2 Cl ({2, 0, 4} - {0, 4, 2} - {1, 5, 3} + {3, 1, 5}) 
+ i e 2 ({0, 1, 3} - {0, 3, 1} - {1, 3, 4} + {3, 1, 4} 
- {0, 2, 3} + {2, 0, 3} - {1, 4, 2} + {1, 2, 4}) . 

Let's observe that, setting the unphysical coefficients e\ and e 2 to zero, T>§ coincides 
with the operator presented in [23]. The coefficient p^l, which can be a priori present 
according to (2.11), actually vanishes. Therefore, the first contribution to the dressing 
phase should be searched at eight loops. We are thus led to consider the eight-loop 
dilatation operator: since it isn't known at present, we now use our procedure to 
constrain it as much as possible. 

The permutation structures contributing to V 8 are listed in Table 1. After imposing 
hermiticity and parity invariance, the eight-loop dilatation operator can be written 
(with real coefficients) as 



V 8 =a { } + V 8 

,even 

+ v 8t0dd + v 8 

,mixed j 



(3.18) 



where 



V s , even = b {0} + c ({0, 2} + {2, 0}) + d {0, 4} 

+ ei ({0, 2, 4} + {4, 2, 0}) + e 2 ({0, 4, 2} + {2, 0, 4}) 

+ *e 3 ({0,4,2}-{2,0,4}) + /{0,6} 

+ g ({0, 2, 6} + {2, 0, 6} + {0, 4, 6} + {0, 6, 4}) 

+ i £l ({0, 2, 6} - {2, 0, 6} - {0, 4, 6} + {0, 6, 4}) (3.19) 
+ Zi({0,2, 4, 6} + {6, 4, 2,0}) 

+ l 2 ({0, 2, 6, 4} + {2, 0, 4, 6} + {0, 6, 4, 2} + {4, 2, 0, 6}) 
+ i e 2 ({0, 2, 6, 4} - {2, 0, 4, 6} + {0, 6, 4, 2} - {4, 2, 0, 6}) 
+ h ({0, 4, 2, 6} + {2, 0, 6, 4}) + m {2, 0, 4, 2} , 



= «i {0, 1} + a 2 ({0, 1, 2} + {2, 1, 0}) + a 3 {0, 3} 
+ a A ({0, 1,3} + {0, 3, 1} + {0, 2, 3} + {2, 0, 3}) 
+ % e ml ({0, 1, 3} - {0, 3, 1} - {0, 2, 3} + {2, 0, 3}) 
+ a 5 ({0, 1, 2, 3} + {2, 3, 0, 1}) + a 6 ({2, 1, 0, 3} + {0, 3, 2, 1}) 
+ % e m2 ({2, 1, 0, 3} - {0, 3, 2, 1}) + a 7 ({0, 1, 4} + {0, 3, 4}) 
+ a 8 ({0, 1, 2, 4} + {4, 1, 2, 0} + {0, 2, 3, 4} + {4, 3, 2, 0}) 
+ i e m3 ({0, 1, 2, 4} - {4, 1, 2, 0} - {0, 2, 3, 4} + {4, 3, 2, 0}) 
+ a 9 ({2, 1, 0, 4} + {4, 1, 0, 2} + {2, 0, 3, 4} + {0, 4, 3, 2}) 
+ i e rn4 ({2, 1, 0, 4} - {4, 1, 0, 2} + {2, 0, 3, 4} - {0, 4, 3, 2}) 
+ a 10 ({0, 1, 3, 4} + {0, 3, 1, 4}) + a u {0, 5} 
+ a 12 ({0,l,5} + {0,4,5}) 
+ £*i 3 ({0, 2, 5} + {2, 0, 5} + {0, 3, 5} + {0, 5, 3}) 
+ % e m5 ({0, 2, 5} - {2, 0, 5} - {0, 3, 5} + {0, 5, 3}) 
+ a u ({0, 1, 3, 5} + {4, 2, 0, 5} + {5, 3, 1, 0} + {0, 2, 4, 5}) 
+ i e m6 ({0, 1, 3, 5} + {4, 2, 0, 5} - {5, 3, 1, 0} - {0, 2, 4, 5}) 
+ a 15 ({0, 3, 1, 5} + {0, 1, 5, 3} + {2, 0, 4, 5} + {0, 4, 2, 5}) 
+ % e rn7 ({0, 3, 1, 5} - {0, 1, 5, 3} + {2, 0, 4, 5} - {0, 4, 2, 5}) 
+ a m ({0, 2, 3, 5} + {2, 0, 5, 3}) + a 17 ({2, 0, 3, 5} + {0, 2, 5, 3}) 
+ *e m8 ({2,0,3,5}-{0,2,5,3}) 
+ ai 8 {0, 3, 6} + ({. . . , a O a + 1, . . .}) . 
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Requiring that the vacuum state is protected we have one condition on the coefficients. 
If we then impose the eight-loop dispersion relation on one-magnon states we obtain 
four more conditions, which in particular fix the coefficient l\ = —5. In order to fix 
other parameters we consider two-impurity states: for sufficiently long states in order 
to avoid wrapping effects, with the help of Mathematica, we explicitly diagonalize the 
dilatation operator and compare its eigenvalues with those obtained from the pertur- 
bative solutions of the asymptotic Bethe equations. For L = 5, 6, 7 (length 10, 12, 14 
states) we find 12 more conditions, while higher length states give no more equations. 
So the point 4 of the procedure described above isn't sufficient to completely fix the 
dilatation operator. Nevertheless, we stress that an important result already follows 
from the equations we could derive from point 4: a non vanishing 6 coefficient (3^1 f° r 
the dressing phase would imply that the ctj coefficients, appearing in T> 8 rnixed , cannot 
be all zero. This proves that T> 8tTnixed ^ 0. In principle, the undetermined coefficients 
could be fixed imposing the commutation of T> with the third charge Q^: this is labo- 
rious and we won't do it, since it isn't necessary for the computations of the dressing 
phase, as we will see in the next chapter. We quote here the part of the dilatation 
operator which acts on even sites only and exchanges the fields inside the operator in 
a maximal way: 

V 8 , mr = + m{2,0,4,2} 

+ h ({0,4, 2, 6} + {2, 0,6, 4}) 

+ ( - 1 + \m + % e 2a + - A f3) ({0, 2, 6, 4} + {0, 6, 4, 2}) (3 21) 

+ ( - 1 + \m - it 2a + ^) ({2, 0, 4, 6} + {4, 2, 0, 6}) 
-5 ({0,2, 4, 6} + {6, 4, 2,0}) 

where the coefficient e 2a corresponds to similarity transformations and depends on the 
particular renormalization scheme used. The dressing phase enters also 2\mr through 
the parameter (5 = fif^. The form of (3.21) is the same of the M = 4 SYM counterpart 
in the four-loop dilatation operator [35] . It is interesting to note that the coefficients ctj 
and the coefficients of h(\) don't appear in (3.21). However, the rest of the dilatation 
operator contains also these parameters: in general, the "mixed" hamiltonian acts 
coupling the two types of magnons; this is a completely new feature of ABJM theory 
which appears at eight loops, making this theory substantially different from the M = 4 
SYM theory at four loops. 

6 We will see in Section 4 that this is, indeed, the case. 
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4. Computation of the Dressing Phase 

In this Section we proceed to the computation of the leading-order coefficient of the 
dressing phase, i.e. the parameter ftfl which appears in the eight-loop dilatation 
operator. We determine it by direct field-theory calculations, in a similar fashion to 
the M = 4 SYM case [35]. Here, we use M = 2 superspace methods for the evaluation 
of Feynman diagrams. We refer to Appendix A for a very brief review of such methods. 

We first rewrite the dilatation operator in the basis of chiral functions introduced 
in Appendix B, which is directly related to the chiral structure of the supergraphs. The 
part of the dilatation operator leading to maximal reshuffling of spins 7 preserves the 
form of (3.21): 

V=mx(2,0,4,2) 

+ h [x(0,4,2,6) + X (2,0,6,4)] 

+ ( - 1 + \m + i e 2a + ^) [ X (0, 2, 6, 4) + X (0, 6, 4, 2)] 

+ ( - 1 + \m - u 2tt + - A f3) [x(2, 0, 4, 6) + X (4, 2, 0, 6)] 
-5(x(0,2,4,6) + X (6,4,2,0)]. 

This part of the dilatation operator is sufficient to compute the parameter fi\ 3 . 2\mr 
has just the same form as the maximal reshuffling hamiltonian of M = 4 SYM [35]. The 
coefficients of the five terms in (4.1) can be computed from the supergraphs contributing 
to x(2,0,4,2), x(0,4,2,6), x(0,2,6,4), x(2,0,4,6), x(0,2,4,6). Only one supergraph 
contributes to each term and it contains only scalar interactions. The coefficient of 
the last term was already fixed by imposing the dispersion relation on one-magnon 
states. The first four terms contain precisely four undetermined parameters, namely 
m, h, t2a and (3: in order to compute such parameters we need to evaluate the super- 
graphs in Figure 1 and isolate the overall UV divergence by subtracting all their UV 
sub divergences. 

First of all, we perform standard D-algebra manipulations (see Appendix A) to 
reduce the supergraphs to ordinary momentum-space integrals I% x . Such integrals are 
displayed in Figure 2. The Feynman rules imply that the color factor is always N 8 , 
which combines with (^r) 8 to give the right power of the 't Hooft coupling, i.e. A 8 
(the (47r) 8 factor will be simplified by the factor present in each integral I 8x )- The 
pole parts of the integrals with the sub divergences subtracted are denoted by I 8x . The 
maximal reshuffling supergraphs, thus, evaluate to: 

7 We give here only the terms which act non-trivially on even sites. 
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<Sa <Sb <Sc <Sd 

Figure 1: Supergraphs corresponding to the maximal reshuffling terms of the eight-loop 
dilatation operator of ABJM theory. From left to right, and together with their reflected 
diagrams, they are associated to each line in (4.1). The horizontal bar represents the operator 
itself, or analogously, the spin chain. 



A 8 (47r) 8 7 8a |i x(2,0,4,2), 




A 8 (47r) 8 J 8b |i X (0,4,2,6) 






A 8 (47r) 8 7 8c |i x(0,2,6,4) 



A 8 (47r) s 7 8d |i x(2,0,4,6) 



A 8 (47r) 8 / 8e |i X (0,2,4,6). 



The subtracted integrals are computed in Appendix D and are given by: 




(4 



7T 



1 



1 



1 , y 



3072 e 4 192 e 3 48 e 2 e 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



'86 



5 5 
+ 



(4tt) 8 V 6144 e 4 768 e 3 384 e 2 32 e 



'8c 



1 1 

+ 



3 5 
+ 



(4tt) 8 V 2048 e 4 128 e 3 Me 2 192 e 
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Figure 2: Momentum-space integrals obtained from supergraphs S a , Si,, S c , S^, S e after 
completion of the D-algebra procedure. 
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For their evaluation, we used the G eg enbauer polynomials x-space technique (GPXT). 
We briefly review this technique in Appendix C. 

We could analytically compute all the integrals in (4.7) apart from the ^ pole 
of h a , denoted by y in the (4.7): in this case GPXT seems to be quite inefficient, 
so we followed a different strategy. We evaluated the integral numerically through 
Mellin-Barnes techniques. Remarkably, guided by the transcendentality principle, it 
was possible to extract the corresponding analytical result from the numerical one by 
means of the PSLQ algorithm [51]. It is given by 



(4.) 



All the details are given in Appendix D. 

We can write now the final results for the supergraphs contributing to maximal 
reshuffling terms of the dilatation operator: 



Sa 


= A 8 


(- 


s b 








= A 8 
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192 


s d 






S e 


= A 8 


5 

16 ; 



5 



1 



32 8 

^X(0,4,2,6), 
X(0,2,6,4), 
^X(2,0,4,6) 



C(3) x(2,0,4,2) 



(4.9) 
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According to the definition (3.3) of the dilatation operator, the corresponding contri- 
butions to the renormalization factor are obtained multiplying these values by -16. At 
this point, we can find the unknown parameters appearing in (4.1). Our results are 8 : 

m = b - - 2C(3) , h = \ , (3 = 4C(3) , e 2a = ^ . (4.10) 

In particular we have computed the value of the leading order coefficient of the dressing 
phase, which is (3^1 = 4£(3). It is the same value found in J\f = 4 SYM for the (3^1 
coefficient. In the Conclusions we comment on this result. 



5. Conclusions 

In this paper we have investigated the dressing phase of ABJM theory at weak cou- 
pling. It appears in the Bethe equations, as well as in the dilatation operator, starting 
at eight loops. A simple procedure based on the Bethe Ansatz allowed us to con- 
struct the asymptotic dilatation operator in the SU{2) x SU(2) sector up to six loops. 
We verified that the dressing phase is not present up to this order. Furthermore, we 
could largely constrain the eight-loop dilatation operator and show that the two SU (2) 
sectors are coupled starting at this order. In this case we focused on the terms corre- 
sponding to maximal interactions: they contain the unknown parameter fi^l, which is 
the leading order coefficient of the dressing phase. Thanks to superspace techniques 
and with the help of the transcendentality principle, we could directly compute its 
value from Feynman supergraphs. A great simplification follows from the fact that we 
considered maximal interactions: we needed to compute a very small number of Feyn- 
man diagrams (in fact, just three are necessary to completely fix which contain 
only scalar interactions. The result we found is P^l = ^C(3) and coincides with the 
value of the leading coefficient of the dressing phase of Af = 4 SYM. In that situation, 
it was argued [17] that the weak-coupling expansion of /? r , s (A) can be extrapolated 
by analytical continuation of the all-loop strong-coupling coefficient functions [52, 16] 

= e~o c&y +s --- \with 

>) (l-(-ir + *CH) fr 1)(s r[(l( 3 + r + n-3)]r[(l( 3 -r + n-l)] 
r ' s 2(-27r)«r(n - 1) 1 A ; r[(±(s + r - n + l)]r[(±(s - r - n + 3)] ' 1 } 



8 The presence of the factor i in the unphysical ti a seems to spoil hermiticity of the dilatation oper- 
ator, see [35] for further comments on this. Moreover, the coefficient of the chiral function %(0, 2, 4, 6) 
turns out to be -5, as it was already known from the previous section: it, thus, constitutes a non trivial 
check of the whole procedure. 
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to negative values of n: 

oo 

PrM = ~Y. c is n) 9 r+s+n ~\ (5.2) 

n=l 

where g — So, the dressing phase coefficients were conjectured to be 

m = -4T- 2e - i] ■ (5.3) 

{—2) 

In particular, cJ> 3 = — 4£(3), so that the leading order coefficient of the SYM dressing 
phase was {3^1 = 4£(3). It was shown that this guess is indeed correct [35]. 

Since in the Bethe Ansatz proposal of [9] the strong coupling limit of the dressing 
phase is the same of the AdS*,/ CFT4 case, we can suppose, by analogy of the M = 4 
SYM case and with the usual replacement g — > h(X), that in ABJM we have 

00 

^ s (A) = -^c(7)/ i (A)^ + - 1 . (5.4) 

n=l 

For r = 2, s = 3, we see that we should obtain (since = = 0) 

& l3 (A) = -h 3 2 c { ~ 3 2) A 6 - (/44; 3 4) + 3/^4; 3 2) ) A 8 H (5.5) 

We observe that, since the first weak coupling coefficient of h(X) is hi = 1, the leading 
order contribution to ABJM dressing phase is just (3\l = 4£(3), which is precisely the 
same value obtained in J\f = 4 SYM. At higher orders, however, the ci~ s n ^ coefficients 
mix with the non trivial coefficients of h(X) yielding dressing phase coefficients which 
differ from the M = 4 SYM ones. We stress that, in our context, this is merely a 
conjecture: anyway, we computed the coefficient by independent field theory tech- 
niques and we found that the above prediction is indeed correct. It would be interesting 
to further test the validity of (5.4) beyond the leading order, through direct computa- 
tions: hopefully, restricting to the case of maximal interactions, the calculations will be 
simplified and have a chance of being performed with standard field theory techniques. 
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A. ABJM Theory in M = 2 superspace 



In this Appendix we briefly review the Af = 2 superspace formulation of ABJM theory. 
This was first given in [53], but we follow the notations used in [21] which are adapted 
from the ones of [54]. ABJM theory has two Af = 2 vector supermultiplets, V and V, 
with V transforming in the adjoint of the first U(N) factor and V in the adjoint of 
the second U(N) factor of the gauge group. In order to extend the supersymmetry to 
Af = 6, the ABJM action also contains two sets of chiral matter superflelds, Z A and 
Wa with A = 1,2. Z A and Wa transform respectively in the (2, 1) and (1, 2) of the 
global SU (2) x SU (2) flavour group. Moreover, they transform in the bifundamental 
representations (N, N) and (N, N) of U(N) x U(N). The gauge fixed ABJM action in 
Af = 2 superspace reads 



s 4 

47T 



| j d 3 x d 4 9 jf * dt Tr (V D a e- tv D a e tv - V D° 'e^ 'D a e tV ) 
+ J d 3 xd i 9^[z A e v Z A e- V + W B e^W B e- v ) 

+ % - J d 3 x d 2 9 e AC £ BD TrZ A W B Z c W D + J d 3 x d 2 9 e AC e BD Tr Z A W B Z C W D 
+ gauge fixing and ghost terms > . 



(A.l) 

The first line contains the Chern-Simons action, the second line contains the kinetic 
term of the matter fields and their coupling with gauge fields, while the third line is 
the superpotential. 

The three-dimensional, Af = 2 superspace spinor covariant derivatives D a , D a sat- 
isfy the algebra 

{Dq, D/3} = {D Q , D/j} = , {D a ,Dp}=p a p. (A.2) 
The metric e A B for the SU(2) flavor indices is given by 

ei2 = 1 , e 12 = 1 , e AB e CD = S A S B - 5 A 5 B ■ (A.3) 

For the integration over the superspace our conventions are J d 2 9 = ^d a d a , J d 2 9 = 
\d a d a and / d 4 ^ = / d 2 #d 2 ^, such that 



J d 3 xd 2 ^ = J d 3 xD 2 | 0=0 - =o , J d*xd 2 9 = J d 3 xD 2 
J d 3 xd 4 9 = J d 3 xD 2 D 2 



\e=e=o 5 

(A.4) 



10=0=0 
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The #-space 5-function is given by 



5\6-6') = {9-9') 2 {9-9') 2 . (A.5) 

We now give the Euclidean Feynman rules (i.e. we have Wick-rotated to e s in the path 
integral) of the theory, relevant for the computations of the diagrams in Figure 1. The 
chiral field propagators are 

(Z B (p)Z A (-p)) = (W B (p)W A (-p)) = - 9 2 ) , (A.6) 



k p 2 

where diagonality in the gauge indices have been suppressed. The vertices are obtained 
by taking the functional derivatives of the Wick rotated action w.r.t. the corresponding 
superfields. When a functional derivatives w.r.t. the (anti)-chiral superfields is taken, 
factors of (D 2 ) D 2 are generated in the vertices. We need only the quartic superpotential 
vertices: 

ie Ac e BD A- [Tr (B^B b _B c -Bd) - Tr(B c -B b _B^)] (A.7) 

ie AC e BD ^ [Yi{BaB h -B c _B±) - Ti{B c _B h -B^)] , (A.8) 

where the color indices are labeled (a, b, c, d) counter-clockwise starting with the leg in 
the upper left corner. Note also that, in a standard way, one of the (D 2 ) D 2 factors has 
been absorbed into the (anti) chiral integration such that the integration measure of 
the (anti) chiral vertex is promoted to the full superspace measure. We have introduced 
matrices B- and Ba, with underlined a = 1, • ■ • , A^ 2 indices that transform in the 
(N,N) of the gauge group U(N) x U(N). 




B. Permutation structures 



The spin-chains which arise in the study of AB JM theory at i loops in the SU(2) x SU(2) 
sector are long-range deformations of the alternating Heisenberg spin chain with next- 
to-nearest neighbor interactions. The interactions among spins are represented by prod- 
ucts of permutations of next-to-neighbor spins, therefore we introduce the permutation 
structures as [19] 

L 

{di, . . . , a n } = 2j P2j +ai) 2j+ai+2 ■ ■ ■ P2j+a n ,2j+a„+2 , (B-l) 
0=1 
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where the permutation operators exchange spins at sites % and j. Indices are 
understood modulo 2L. If n = 0, we denote { } = L. The range 72. of a permutation 
structure is 

1Z = max(ai, . . . , a n ) — min(ai, . . . , a n ) + 3 . (B.2) 

The integer n is called length of the permutation structure and, in terms of the associ- 
ated Feynman diagrams, it coincides with the number of chiral and antichiral vertices. 
The permutation structures satisfy the following relations: 

{... ,a,a, ...} = {...,...} ; 

{... ,a,b, ...} = {... ,b,a, ...} , if \b — a\ ^ 2 ; 

{a,b, ...} = {a + 2m,b + 2m, ...}, m = 0,1,2,... (B.3) 
{. . . , a, a + 2, a, . . .} = {. . . , a, a + 2, . . .} + {. . . , a + 2, a, . . .} - {. . . , a, . . .} 

-{...,a + 2,. ..} + {.. .,...}. 

Under hermitian conjugation and parity, they transform as 

{ai, . . . ,a n y = {a n , ...,a 1 }; 
V{a 1 ,...,a n }V = {-ai, . . . , -a n } . 

The permutation structures represent a convenient basis for the dilatation operator. 
Nevertheless, when dealing with Feynman diagrams, it is useful to use a different basis, 
namely the basis of chiral functions, defined in terms of permutation structures as 
[55, 56, 50] 

*()={}, 
x(a)={a}-{}, 

x (a,b)={a,b}-{a}~{b} + {}, 

X (a, b, c) = {a, 6, c} - {a, b} - {a, c} - {b, c] + {a} + {b} + {c} - { } , (B.5) 

x(a, 6, c, d) = {a, 6, c, d} — {a, 6, c} — {a, 6, d} — {a, c, rf} — {6, c, d} 

+ {a, 6} + {a, c} + {a, rf} + {b, c} + {6, d} + {c, c?} 

-W-W-W-{4 + {}- 

These functions are directly related to the chiral structure of a supergraph and precisely 
describe the flavor flow of fields inside (2.1) under the action of the interaction. All 
diagrams contributing to the same chiral function share the same chiral structure. For 
the maximal interactions encountered in this work, no vector interactions can be present 
and the chiral structure alone determines the unique diagram which contributes. The 
chiral structure of the function x(°i> • • • , a n) contains n chiral vertices and n antichiral 
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vertices, connected by (ZZ) or (VFW^) propagators. Each pair of chiral and antichiral 
vertices (plus the propagator connecting them) describes an elementary permutation 
of next-to-nearest neighbors. 



C. The Gegenbauer Polynomial x- Space Technique 

In this Appendix we review the Gegenbauer Polynomial x-Space Technique (GPXT), 
which has been applied to find the UV divergent part of the Feynman integrals needed 
in the present work. The technique was introduced in [57] and developed in [58, 59]. 
See also [60, 50] for useful reviews. All the integrals are computed using dimensional 
regularization in Euclidean space of dimension 

£> = 2(A + 1), (C.l) 

where 

A = --e (C.2) 
in order to get three-dimensional space in the e — > limit. 

In our computation, a generic Wick-rotated £-loop integral in momentum space has 
always the form 

1 f d D k x ■ ■ ■ d D k e 



Ie -J^r J Hi---np ' (C ' 3) 

where P is the number of propagators IT, which, in general, depend on the loop mo- 
menta ki, . . . , ke and the external momenta pi, . . . ,p e . We always consider massless 
propagators. 

A propagator with weight a in momentum space is Fourier-transformed to coordi- 
nate space according to the following formula 

j_ r(A + l-a) f d D xe^ 

k 2a r(a)7r A+1 J x 2 ^ 1 -^ ' 1 ' ' 

where T(z) is the Euler gamma function. The weight a is a generic complex number. 
We stress the presence of the unconventional factor 2 in the exponential. According to 
this definition, we have 

' d D ke 2tkx = n 2 ^5(x). (C.5) 



According to GPXT technique, the computations are made directly in coordinate space 
rather than in momentum space. The technique is grounded on the observation that, 
in x-space, the scalar propagator always depends on the difference of two points, 

A(xi, Xj ) = , (C.6) 
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and can thus be expanded in terms of the Gegenbauer polynomials, which form an 
orthogonal set on the unit sphere in R D . For the moment, D is an arbitrary integer 
dimension and A = D/2 — 1. The analytic continuation to complex D will be done in 
a second step. 

The Gegenbauer polynomials C A are defined in terms of a generating function, 

oo 

5>C*(*)f B , (C7) 



;i - 2xt + t 2 y 



n=0 



where x G [—1, 1]. We refer to the quantity A as the weight of the polynomial, while 
n is its index. The Gegenbauer polynomials are orthogonal with respect to the weight 
function (1 — x 2 ) x ~ 1 / 2 : 



7r2 1 - 2A r(n + 2A) 
, n\(n + X)T(\y 



dx (1 - xy-C^x)C x m (x) = - ■ \— — 5 nm . (C 



The following particular values are often needed: 

T(n + 2a) 



c a n {i) =- 



n!r(2a) ' 



(C.9) 



Let's now turn to the description of the technique. We will consider here only 
diagrams with a single external momentum p, entering and leaving the graph at points 
x in and x out in x-space. The generic £-loop integral in (C.3) becomes, in x-space, 

TVX'l-P P 2ip(x out -x in ) 

' 1 d D Xl ---d D x P J- — (C.IO) 



1 (¥tt p ) x+1 J Ax 

where Aj are the propagators in coordinate space. 

It is now convenient to move to spherical coordinates in D dimensions: to this 
purpose we define 

r = x\ x = ^=. (Gil) 



So r is the (squared) radial coordinate and x is the unit vector pointing in the same 
direction as x. The integration measure changes to 

d D x = ^S D ^r x drdx, (G12) 

where 

2tt a+1 

So-, = f ^ (C.13, 
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is the surface of the unit sphere in R . The integral (CIO), therefore, transforms to 

If = Nx y p) f dri ■ ■ ■ drp^dx! ■ ■ ■ dx P ^r x ■ ■ ■ r x P „ f e *p(*°ut-*m) 
J Ai ■ ■ ■ A P 

with the normalization factor 

W = |r|' (ai5) 

At this point we can expand the propagators in terms of the Gegenbauer polynomials: 
from (C.7) we have 

= J^W = W, • (C16) 

where we have introduced the notation 

rriij = min(r i? r,-) , Mj ;j = max(r;, rj) . (C17) 

If we are interested only in the UV divergent part of the loop integral, as in all the cases 
in this work, a great simplification occurs: in x-space, UV divergences appear when 
all coordinates are small, so we can approximate the exponential factor with unit, ne- 
glecting it. Since it is equivalent to set to zero the external momentum p, dropping the 
exponential will introduce IR divergences, in the region where some coordinates are 
large, which mix to UV ones altering the final result. We thus regulate them by intro- 
ducing an infrared cutoff R in the radial integrations. At the end of the computation, 
when all the sub divergences have been subtracted, the principal part of the integral 
must be independent of the regulator R, while the finite part in general depends on R 
and should be discarded from the result. 

The expansion of the propagators in series of Gegenbauer polynomials allows us 
to separate the integral in radial and angular parts and introduces as many infinite 
sums as the number of the propagators themselves. However, we can minimize the 
number of series in this way: thanks to translational invariance of the integral (CIO), 
we can choose one of the vertices of the diagram as the origin of D-dimensional space. 
We call this vertex the root vertex. All the propagators directly connected to the root 
vertex, which we call root propagators, are simply given by l/r^, where is the radial 
coordinate of the vertex which connects the propagator to the root vertex, and so they 
don't produce any series expansion in the Gegenbauer polynomials. Obviously, the 
best choice of the root vertex is usually such that the number of propagators attached 
to it is maximized. Therefore, in most cases, the root vertex will coincide with the 
composite operator. 
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Once we have chosen the root vertex and we have expanded the propagators de- 
pending on differences of coordinates in terms of Gegenbauer polynomials, the angular 
and radial integrations are performed separately. 

The angular integration can be performed by repeated use of the orthogonality 
relation (C.8) of the Gegenbauer polynomials, which can be rewritten in terms of the 
angular variables as 



The angular integration is normalized as f dx — 1. In particular, since Cq(x) = 1, we 
have 



Kronecker deltas can be used to decrease the number of summations. Angular 
integration is usually the hardest part in multiloop computations and can be extremely 
simplified if the root vertex is chosen so as to minimize the angular loop number and 
the number of infinite summations. In fact, GPXT is at its best when the number of 
infinite summations can be reduced at most to one. 

Let's now turn to the radial integration. Having dropped the exponential and 
introduced the infrared cutoff, radial integrands consist of simple powers. The only 
difficulty is that, because of the presence of the min and max functions, the domain of 
integration (which is an hypercube of length R) has to be split into (P—£)\ subdomains, 
defined by the different orderings of the radial variables. This number can be large. 
Of course this is not a problem if the procedure is automated with the help of a 
computer 9 . Anyway, it is useful to find all the possible symmetries of the integrand in 
order to reduce the independent domains of integrations. 

At this point, when the angular and radial integrations have been performed, the 
next point is to promote the dimension D, or equivalently A, to a complex parame- 
ter through the formula (C.2) and then perform the Laurent expansion of the result 
around e — 0. If multiple poles are present, we proceed to the subtraction of subdi- 
vergences. These must be computed within the same renormalization scheme as the 
original integral, i.e. using GPXT and introducing the same cutoff procedure for in- 
frared regularization. 

The last step is to perform the summations that possibly survived after the angular 
integrations. As stated before, finding analytical results can be very hard, especially 
when multiple series are present, and sometimes only a numerical analysis is possible. 

9 In particular, we used the Mathematica routine for radial integrals described in [60]. 




(C.18) 




(C.19) 
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D. Integrals 



In the following, we list the results of the integrals relevant for the computation of 
maximal reshuffling diagrams of AB JM theory in the SU(2) x SU (2) sector (see Figure 
1). They are computed via GPXT as described in Appendix C. We denote by Ij 
the value of the integral. We give here the e expansions up to the order needed for 
the computation of the UV divergences of the diagrams in Figure 1. We denote by R 
the infrared regulator. For convenience, we leave here the normalization factor (C.15) 
unexpanded: this will be easily reintroduced in the subtraction of subdivergences. 



iV A (2, 3) ( + \ogR + e log 2 R+\e 2 log 3 R 



(D.l) 
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1 1 + log R 



2e 



- l + 21og J R + log 2 J R 



= iV A (6,9) 
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= iV A (4,6)(^ + -^ + 21og 2 i?+^log 3 /?), (D.3) 
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(D.9) 



(D.10) 



(D.ll) 
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where, in I 8a , x is reduced to a multiple series which we don't show here (see below 
for a more effective method to compute the Ig a integral). The subtracted integrals 
are denoted by Ij and are graphically represented by a box around the corresponding 
integral. We conveniently factor out l/(47r) £ from the results. The pole parts of such 
integrals are given by: 



1 1 

(4tt) 2 47 
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(D.22) 
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Regarding I$ a , GPXT allowed us to analytically compute the higher order poles in e 
and to reduce the first order pole y to a multiple series. Such a series is not easy to sum 
and we choose to compute the integral Ig a through Mellin-Barnes technique [61, 62] in 
order to have a reliable numerical estimate of its pole Is a \i- 

First of all, we contract the bubbles to reduce to the evaluation of a four-loop 
integral: 





G(l,l) 3 G(l,l/2 + e) 



where the "G-functions" are defined by 

T(A + 1 - a)r(A + 1 - P)T(a + (3 - X - 1] 



(D.25) 
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(47r) A + 1 r(a)r(/3)r(2A + 2 - a - f3) 
Since the G-functions in (D.25) present a simple pole in e, the four-loop integral 
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has to be computed up to order zero in e. Let's introduce the corresponding four-loop 
master integral with generic powers of the propagators: 



J 46 (ai, . . . ,a 8 ) 



1 




d D k x d D k 2 d D hd D k A 



1 



i(h - h- P ) 2 ]^[{k 2 - fa- P )*]°°[(ki - MTIOs - h) 2 ]^ ' 

(D.28) 

The Mellin-Barnes representation of this integral and its analytical continuation in e 
can be obtained with the help of the Mathematica packages AMBRE [63] and MB [64] 
respectively. 

The MB package also contains a routine suitable for numerical integration of MB 
representations. Besides the built-in Mathematica function NIntegrate, it uses the 
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CUBA library [65] of numerical integration routines and the CERN libraries [66] for the 
implementation of gamma and psi functions, in order to prepare Fortran programs, 
which are more efficient, in terms of computational time, for high dimension MB inte- 
grals. 

A MB representation for the integral J 4b is 

dz 6 r(A + 1 - a 17 - zi)r(-zi)r(A + 1 - a 15 - z 2 ) 



X 



X 



X 



X 



2ni r(a 1 )r(a 3 )r(a 5 )r(« 6 ) 

r(-* 2 )r(ai + z 12 )T(X +l-a 26 + Zl - z 3 )T(2X + 2 - a l257 - z 2A )T(-z 4 ) 
T(2A + 2 - a 157 )r(a 7 )r(a 2 - Zl )F(3X + 3 - a 12567 - z 2 ) 

T(a 2 - z 1 + z u )T(-2X -2 + a 12567 + z 2M )T(3\ + 3 - ai 2 35678 ~ ^2345) 
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T(-3A - 3 + ai 2356 78 + z 5 q) 

(D.29) 

where we have denoted a^k... = + cxj + + • • • and similarly for z^. 
We have to compute 

hb = Ml/2 + £ , 1, 1, 1/2 + e,2e, 1,1,1). (D.30) 

Putting this integral into the MB package we can make the ^-expansion and find the 
following numerical result: 

2.116213934935895 1(T 8 (8.26225 ± 0.00003) 1(T 7 

1 46 — o 1 5 



(D.31) 



(0.0000177432 ± 0.0000000003) , 
+ - - + (0.0002705914 ± 0.0000000003) . 



We then insert this result in (D.25), make the e-expansion and subtract the subdiver- 
gences, to obtain the following numerical result: 

- 5.2348 . . . 10~ 13 8.37567 . . . 10~ 12 3.35028 . . . lO" 11 9.63613 . . . 10~ 12 
ha -, V 



e 4 e 3 e 2 e 



(D.32) 
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Comparing this result with (D.20) we see that the poles of second, third and fourth 
order perfectly agree. We recall that these poles were computed analytically via GPXT. 
In addition, we now have the numerical result for the first order pole: 

I 8a \i/e = (-0.0059921 ± 0.0000008) . (D.33) 

We now want to extract the analytical result from this number. In section 2 we argued 
that the parameter /3 of (4.1) is purely transcendental and that it can be a rational 
combination of just two transcendental constants, namely n 3 and C(3). Since the only 
transcendental contribution to (3 can come from the ^ pole of I$ a , such pole should be a 
rational combination of the constants 1, 7r 3 , £(3). In order to extract the coefficients of 
this generic linear combination from the numerical value (D.33) we use the Mathematica 
implementation [67] of the PSLQ algorithm 10 [51] (see also [68]). This is a powerful 
integer relation detection algorithm: given a vector x = (xi, . . . , x n ) of real or complex 
numbers and a precision 10~ p , PSLQ looks for a vector a = (ai, . . . ,a n ), if exists, of 
integers a«, not all zero, such that 

aiXi H h a n x n = . (D.34) 

In our case, the vector x contains the numerical value in Ig a |i/ e and the transcendental 
constants it should be fitted with: 

x = (-0.0059921, 1,tt 3 ,C(3)). (D.35) 

Application of PSLQ algorithm with precision 10~ 7 , which is the same within which 
the numerical value of I$ a |i/ e has been computed, gives the solution 

a = (-32,-5,0,4) . (D.36) 

It means that the exact value of the ^ pole of the I$ a integral should be 

W= (S5i(-| + 5«3)). (D.37) 
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